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I. INTRODUCTION 


Choptuik’s discovery of critical phenomena in the gravitational collapse of a scalar held 
m] sparked a snrge of interest in gravitational collapse jnst at the threshold of black hole 
formation. The discovery of critical behavior in several other matter models qnickly followed 
M- Despite the fact that the evolution equations are very complex and highly non-linear, 
the dynamics of the near-critical held evolution is relatively simple and, in some important 
aspects, universal. The critical solution, which depends on the matter model only, serves 
as an intermediate attractor in the phase space of solutions, and often has an additional 
peculiar symmetry called self-similarity. The mass of the black hole produced in supercritical 
evolution scales as a power law 


Mbh(p) oc (1) 

with parameter p describing initial data, and mass-scaling exponent (3 is dependent only 
on the matter model, but not on the initial data family. An interesting consequence of 
mass scaling which has direct bearing on the cosmic censorship conjecture is the fact that 
arbitrarily small black holes can be produced in near-critical collapse, with the critical 
solution exhibiting a curvature singularity and no event horizon. 

The explanation of the universality of the critical behavior lies in perturbation analysis 
and renormalization group ideas Mi- It turns out that critical solutions generally have 
only one unstable perturbation mode, making them the most important solutions for under¬ 
standing the dynamics of held evolution, after the stable ones (hat space and Schwarzschild 
or Kerr-Newman black hole). As the near-critical held conhguration evolves, all its per¬ 
turbation modes decay, losing information about the initial data and bringing the solution 
closer to critical, except the one growing mode which will eventually drive the solution to 
black hole formation or dispersal, depending on its content in the initial data. Thus the 
critical solution acts as an intermediate attractor (of codimension one) in the phase space 
of held conhgurations. Finding the eigenvalue of the growing perturbation mode allows one 
to calculate important parameters of the critical evolution, the mass-scaling exponent in 
particular. 

An important question is how generic the critical behavior is with respect to initial data, 
or, in phase space language, how big is the basin of attraction of the critical solution. So far 
most of the work on critical gravitational collapse, numerical or analytic, has been restricted 
to the case of spherical symmetry, simply because of the enormous difficulties in treating 
fully general non-symmetric solutions of Einstein equations. A natural concern is whether 
the critical phenomena observed so far are limited to spherical symmetry, and whether 
the evolution of non-spherical data will lead to the same results. The numerical study of 
Abrahams and Evans on axisymmetric gravitational wave collapse |0] and recent numerical 
perturbation calculations by Gundlach i,0 give numerical evidence for the claim that 
critical phenomena are not restricted to spherical symmetry, and that the critical solutions 
are indeed attractors in the full phase space. In this paper we search for analytical evidence 
to support that claim. 

One of the few known closed form solutions related to critical phenomena is the Roberts 
solution, originally constructed as a counterexample to the cosmic censorship conjecture 
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0, and later rediscovered in the context of critical gravitational collapse [^,0. It is 
a continnously self-similar solntion of a spherically symmetric gravitational collapse of a 
minimally conpled massless scalar held. While it is not a proper critical solntion, as it has 

it is still a good (and simple) toy model of the critical 
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more than one growing mode 
collapse of the scalar held. 

This paper considers fnlly general pertnrbations of the Roberts solntion in a gange- 
invariant formalism. Dne to the symmetries of the backgronnd, the linear pertnrbation 
eqnations deconple and the variables separate, so an exact analytical treatment is possible. 
We hnd that there are no growing pertnrbation modes apart from spherically symmetric 
ones described earlier . So all the non-sphericity of the initial data decays in the collapse 
of the scalar held, and only the spherically symmetric part will play a role in the critical 
behavior. 

To onr knowledge, this is the hrst paper to obtain analytical resnlts on non-spherical 
critical collapse. 


II. THE ROBERTS SOLUTION 


The spacetime we will nse as a backgronnd in onr calcnlations is a continnonsly self¬ 
similar spherically symmetric solntion of the gravitational collapse of a massless scalar held 
(the Roberts solntion). The Einstein-scalar held eqnations 

Rfiu = ( 2 ) 

□0 = 0 (3) 

can be solved analytically in spherical symmetry by imposing continnons self-similarity on 
the solntion, i.e. by assnming that there exists a vector held ^ snch that 


£^(f) 0 , 


( 4 ) 


where £ denotes Lie derivative. Self-similar solntions form a one-parameter family, which 
is most easily derived in nnll coordinates |12, 1^,15 
metric 


The critical solntion is given by the 


ds^ = —2 du dv -|- 


( 5 ) 


where 


r = 




— uv, 


0 = 



V 

u. 


( 6 ) 


The global strnctnre of the critical spacetime is shown in Fig. |^. The inhux of the scalar 
held is tnrned on at the advanced time n = 0, so that the spacetime is Minkowskian to the 
past of this snrface. The initial conditions for the held eqnations (U) and (||) are specihed 
there by the continnity of the solntion. 

It is instrnctive to rewrite Roberts solntion in new coordinates so that the self-similarity 
becomes apparent. For this pnrpose we introdnce scaling coordinates 
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s = — ln(— m), 


( 7 ) 



V 

ui 


with the inverse transformation 


u = 


V = 


!)• 


( 8 ) 


The signs are chosen to make the arguments of the logarithm positive in the region of interest 
(n > 0, M < 0), where the held evolution occurs. In these coordinates the metric (H) becomes 


[(1 - e-^^)ds‘^ - 2dsdx 
and the critical solution is simply 




( 9 ) 


r = e 


= X. 


( 10 ) 


Observe that the scalar held 0 does not depend on the scale variable s at all, and the only 
dependence of the metric coefficients on the scale is through the conformal factor This 
is a direct expression of the geometric requirement in scaling coordinates; the homothetic 
Killing vector ^ is simply — 


III. GAUGE-INVARIANT PERTURBATIONS 


To avoid complicated gauge issues of fully general perturbations, we will use the gauge- 
invariant formalism developed by Gerlach and Sengupta jlblJT^. This formalism describes 
perturbations around a general spherically symmetric background 


g^^dx^dx'^ = QAsdx^dx^ + r‘^jabdx°'dx^ 


( 11 ) 


which in our case we take to be the Roberts solution (^. Here and later capital Latin indices 
take values {0,1}, and lower-case Latin indices run over angular coordinates, qab and r are 
dehned on a spacetime two-manifold, while 'jab is the metric of the unit two-sphere. 

Because the background spacetime is spherically symmetric, perturbations around it 
can be decomposed in spherical harmonics. Scalar spherical harmonics Yim{0^ip) have even 
parity under spatial inversion, while vector spherical harmonics Sima{d,(p) = eJ^Yim,b have 
odd parity. We will only concern ourselves with even-parity perturbations here, since odd- 
parity perturbations can not couple to scalar held perturbations. We will focus on non- 
spherical perturbation modes (/>!), as the spherically symmetric case {I = 0) was studied 
earlier 0- For clarity, angular indices l,m and the summation over all harmonics will be 
suppressed from now on. The most general even-parity metric perturbation is 


Sgaudx^dx'^ = hAsYdx^dx 


B 


+hAY^b{dx^dx^ + dx^dx^) 

+ GY:Adx‘‘dx\ 


( 12 ) 


and the scalar held perturbation is 
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5(t) = FY. 


( 13 ) 


As you can see, metric perturbations are described by a two-tensor Hab, a two-vector Ha, 
and two two-scalars K and G; the scalar field perturbation is described by a two-scalar F. 
However, these perturbation amplitudes do not have direct physical meaning, as they change 
under the (even-parity) gauge transformation induced by the infinitesimal vector field 

^,dx^ = UYdx^ + ^Yadx\ (14) 

One can construct two gauge-invariant quantities from the metric perturbations 


kAB — hAB — 2p(A|B), 

k = K — 2 v^pa-, (15) 

and one from the scalar field perturbation 

f = F + 4^’^PA, (16) 


where 


Va = 




Pa = hA- 



(17) 


Only gauge-invariant quantities have physical meaning in the perturbation problem. All 
physics of the problem, including the equations of motion and boundary conditions, should 
be written in terms of these gauge-invariant quantities. Once gauge-invariant quantities 
have been identified, one is free to convert between gauge-invariant perturbation amplitudes 
and their values in whatever gauge choice one desires. 

We will work in longitudinal gauge (h^ = G = 0), which is particularly convenient since 
perturbation amplitudes in it are just equal to the corresponding gauge-invariant quantities. 
The above condition fixes the gauge uniquely for non-spherical modes. (There is some gauge 
freedom left over for the I = 0 mode, but remember that we are only concerned with higher 
I modes.) Expressions for the components of the linear perturbation equations 


(18) 

5(n0) = 0 (19) 

for a fully general perturbation in longitudinal gauge are collected in Appendix By 
inspection of the Op component of the equations, it is clear that the equations of motion 
require that h^v = 0 for / > 1. With the change of notation huu = U and h^v = V, the 
remaining equations of motion for non-spherical modes are 

4(u^ — uv)F^yu — uU^y — uK^u + vK^y -|- vV^u — 2uF^y + 2{2u — v)F^y -|- 2l{l + 1)F = 0, (20a) 


- 2(m^ - uv)K^uu + uU^u + (2n - v){U^y - 2K^y) - 4nF„ -F l{l l)U = 0, (20b) 


(n UV^iUyy “1“ 2K yy Y Yyy^ Y UU y Y UF 


{2u - v){K^y Y Vu) + ‘^uF^u - “^vF^y = 0, 

(20c) 
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— 2{u^ — uv)K^vv + 2uK^y — uVu — (2m — v)Vv + 4mF^ + l(l + 1)V = 0, (20d) 

2(m^ — uv)K^^u — uU^y — 2uK^u + (2m — v){2K^^ + Vy) — ‘^K + /(/ + 1)K + 2V = 0, (20e) 

(m^ — uv){U^y + K^u) + 2,vF = 0, (20f) 


(m - v){V^u + K^y) - 2F = 0. 


(20g) 


Equation (|20aj) comes from the scalar wave equation, and equations ( P0b| - P0^) are the mm, mm, 
MM, 99, u9, and v9 components of the Einstein equations, correspondingly. As usual with 
a scalar held, the system (pOf ) has one redundant equation, so equation (|20c|) is satished 
automatically by virtue of other equations. Equations 


and (|20g]) are constraints, and 
the remaining four equations are dynamic equations for four perturbation amplitudes U, V, 
K, and F. 

Boundary conditions for the system (^) are specihed at m = 0 and the spatial inhnity. 
Continuity of matching with hat spacetime at the hypersurface m = 0 requires the vanishing 
of the perturbations there. We also require well-behavedness of the perturbations at I~ and 
X’*', so that the perturbation expansion holds. Thus, the boundary conditions are 


U = l/ = A = X = 0atM = 0, 

U, V, K, F are bounded at m = — cxd and v = -l-oo. 


( 21 ) 


Equations ( pOD together with boundary conditions ( pT]) constitute our eigenvalue problem. 


IV. DECOUPLING OF PERTURBATION EQUATIONS 


It is possible to decouple the dynamic equations (20a- |20j) by combining them with 
the constraints ([201|) and (|2Ug| ), and their hrst derivatives. After somewhat cumbersome 
algebraic manipulations, which we will not show here, the system of linear perturbation 
equations o can be rewritten as 


2vF 


2(u^ — uv)F^yy — uF^u + (2m — v)F^y H- h l{l + 1)F = 0, 


M — M 


(22a) 


2(m^ — uv)U^yu + uU^u + 3(2m — v)U^y -|- l{l -|- 1)17 = 0, 


2(m^ — uv)Vyu — 3uVu — (2m — m)17.u -|- -|- 1)17 = 0, 


2(m^ — uv)K^yu — uK^u + (2m — v)K^y — 2K + + 1)K = —217 — 


AuF 


u — V 


(22b) 

(22c) 

(22d) 


2vF 

UU y UK y H-= 0, 

M — M 


(22e) 
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2F 

Vu + K,, -= 0. (22f) 

u — V 

This decoupled system of partial differential equations can be further simplihed by exploiting 
continuous self-similarity of the background to separate spatial and scale variables. With 
this intent, we rewrite equations in terms of the scaling coordinates 

1 (1 _ + F,, + F, - 2(1 - e-2")F + + 1)F = 0, (23a) 

1 (1 _ - 217,. - t/,. + l{l + 1)17 = 0, (23b) 

— (1 — e ^*)h(xx -|- 17.S + 217. -|- 317s -|- 1(1 -|- 1)17 = 0, (23c) 

i (1 - e-2")F.. + F.S + Fs -2K + 1(1 + 1)F = -21" - de'^'F, (23d) 

f/,. - (1 - e2^)F. + 2e2^Fs - 4(1 - e^^)F = 0, (23e) 


F. - (1 - + 26^^17, + 4F = 0. 


(23f) 


We decompose the perturbation amplitudes into modes that grow exponentially with the 
scale s (which amounts to doing Laplace transform on them) 

F{x,s)=Y.F,{x) e"^ 

K 

U{x,s)=J2U^{x)e^\ 

K 

V{x,s)=Y.V^{x)e>^\ 

K 

K{x,s)=J2K^{x)e^\ (24) 


The summation runs over the perturbation mode eigenvalues k, which could, in general, be 
complex. Modes with Re k > 0 grow and are relevant for critical behavior, while modes 
with Re/t < 0 decay and are irrelevant. The growing perturbation mode amplitudes vanish 
at s = —oo, so the boundary condition at I~ is satished automatically. For clarity, the 
perturbation mode subscript k and the explicit summation over all modes will be suppressed 
from now on, so henceforth F, U, V, and F will mean F^, Uf^, V^, and K,^ for the mode 
with eigenvalue k. 

The decomposition (|^) converts the system of partial differential equations (|^) into a 
system of ordinary differential equations, which is much easier to analyze: 


- (1 - e-^^)F" + kF' + kF- 2(1 - e-^^)F + l{l + 1)F = 0, 


(25a) 
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^ (1 - e-‘^^)U" + (k - 2)U' -kU + l{l + l)U = 0, (25b) 

^ (1 - e-‘^^)V" + {k + 2)V' + 3kV + l{l + l)V = 0, (25c) 

^ (1 - e-‘^^)K" + kK' + {k- 2)K + /(/ + l)K = -21^ - 4e"^^F, (25d) 

U' -{I- e^^)K' + 2Ke^^K - 4(1 - e^^)F = 0, (25e) 

/s:'-(l-e^^)\/' + 2Ke2^V^ + 4F = 0. (25f) 

The prime denotes a derivative with respect to spatial variable x. These equations can be 
converted into standard algebraic form by the change of variable 

y = e^^, x=]^\ny, (26) 

so that the system (|25|) becomes 

2/(1 - 2 /)$ + [3 - (k + 3)2/]<i> - [3k/2 + l{l + l)/2]4> = 0, (27a) 

2/(1 -y)U+[l-{K- l)y]U - [-«:/2 + /(/ + l)/2][/ = 0, (27b) 

2/(1 - y)V + [1 - (k + 3)y]V - [3 k/2 + + l)/2]V = 0, (27c) 

2/(1 - y)k + [1 - (k + l)y]k - [k/2 -! + /(/ + l)/2]K = 2^ + V, (27d) 

U + {y — l)k + kK + 22/4* — 24> = 0, (27e) 

iir +(2/-1)14 + ^17 +24> = 0. (27f) 


The dot denotes a derivative with respect to 2/, and we redehned the scalar held perturba¬ 
tion amplitude as F = y^ to cast the equations into standard table form. The boundary 
conditions (^) are 

U = V = K = ^ = 0aty = l, 

U, V, K, 2/4) are bounded at y = -|-cx). (28) 

Imposed on the decoupled system of ordinary differential equations (^Tj), these boundary 
conditions give an eigenvalue problem for the perturbation spectrum k. 
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V. PERTURBATION SPECTRUM 


In the previous section we formulated an eigenvalue problem for the spectrum of non- 
spherical perturbations of the critical Roberts solution. We now proceed to solve it. Observe 
that equations (|27a| - ^7dD governing the dynamics of the perturbations are hypergeometric 
equations of the form 


y{l — y)X + [c — (a + 6 + l)y]X — abX = 0. 


(29) 


Equation ( p7d|) is not homogeneous, but we will deal with that shortly. The hypergeometric 
equation coefficients are different for equations describing the perturbations d*, U, V, and 
K] they are summarized in the table below. 



c 

a + b 

ab 

$ 

3 

K + 2 

+ ^l{l + 1) 

u 

1 

K — 2 

— + ^l{l + 1) 

V 

1 

K + 2 

|/t) + ^/(/ + 1) 

K 

1 

K 

+ ^/(/ + 1) 


Hypergeometric equations have been extensively studied; for complete description of their 

. Hypergeometric equation (^) has three singular points at 
y = 0 , 1 , oo, and its general solution is a linear combination of any two different solutions 
from the set 


properties see, for example. 


Xi = IF {a, b] a + b + 1 — c; 1 — y), 

X 2 = {1 — yY~°‘~'^!F{c — a, c — b] c + 1 — a — b-,1 — y), 

^3 = a + l-c;a + l-6; y~Y, 

-^4 = + 1 - c, 6; 6 + 1 - a; y~Y, (31) 


where lF{a,b-, c;y) is the hypergeometric function, which is regular at y = 0 and has 
jF(a, 6;c;0) = 1. Any three of the functions (pi|) are linearly dependent with constant 
coefficients. In particular. 


. _ r(c+l-a-6)r(6-a) , 

^-r(i-o)r(c-a)—' 

r(c + 1 - a - b)T{a - b) ^ 

r(i - (,)r(c -1) ' 


(32) 


The functions Xi, X 2 are appropriate for discussing the behavior of solution near y = 1, 
while X 3 , X 4 give the behavior at inhnity 


Ai = 1, X 2 = (1 - yY “ ^ near y = 1, 

-^3 = (-2/)"“, -^4 = (-2/)“^ near y = 00 . (33) 
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As we said before, imposing the boundary conditions (^) on solutions of equations (^7[) 
leads to a perturbation spectrum. We will now investigate what restrictions the boundary 
conditions place on the hypergeometric equation coefficients. The vanishing of perturbation 
amplitudes aiy = 1 rules out Xi as a component of the solution and requires that Re(c —a — 
6 ) > 0 to make X 2 go to zero. The solution X 2 has non-zero content of both X 3 and X 4 by 
virtue of (l3^) , hence for it to be bounded at inhnity, both Rea and Re 6 must be positive to 
guarantee convergence of X 3 and X 4 . So, unless there is degeneracy, the boundary conditions 
translate to the following conditions on the hypergeometric equation coefficients: 


Re(c — a — 6 ) > 0, 
Rea, Re 6 > 0. 


(34a) 

(34b) 


We are now ready to take on system (|^) . Take a look at equation (|27d) for V. Condition 
for it is ReK < —1, i.e. there are no growing V modes! With the amplitude of relevant 


(36) 


V perturbation modes being zero, the constraints (|27e|) and (P71|) become 


K = -, <h = —, 

K 2k 


and right hand side of equation (|27d|) can be absorbed by the left hand side, making the 
equation for K homogeneous (with c = 2). Indeed equations (p7d|) and (|27a|) for K and d* 
are just derivatives of equation ( P7b|) for U 

(36) 


2/(1 - y)U +[!-(«- l)y]U - [-«/2 + l{l + 1)/2]U = 0, 
which is the homogeneous hypergeometric equation with coefficients 

c=l, a + b = n — 2, ab = 


(37) 


Imposing the boundary condition ai y = 1 for the solution of the above equation and its 
derivatives, which behave like 


f/ cx (1 — y)^~'^ 
K (X {1 — yY~'^ 
$ cx (1 — yy~'^ 


near y = I, 


(38) 


produces restriction on the non-spherical mode eigenvalue 

Re K < 1, (39) 

which is the strongest of restrictions ( 34a ) for equations for U, K, and <h. But then 

Re a -|- Re 6 = Re K — 2 < —1, (40) 


and hence Re a and Re b can not be both positive, and so the boundary condition at inhnity 
can not be satished. A more careful investigation of degenerate cases of relation (^) shows 
that the contradiction between boundary conditions at y = 1 and inhnity still persists if 
R = 0. It can only be resolved by the trivial solution U = K = ^ = 0. Thus we have shown 
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that there are no growing non-spherical perturbation modes around the critical Roberts 
solution. 

In fact, an even stronger statement is true. The contradiction between boundary condi¬ 
tions sX y = 1 and inhnity can not be resolved by a non-trivial solution so long as R = 0, 
i.e. so long as Re/t > —1. Hence non-trivial non-spherical perturbation modes of critical 
Roberts solution must decay faster than e~^. 


VI. CONCLUSION 


In this paper we used the gauge-invariant perturbation formalism to explore the crit¬ 
ical behavior in the gravitational collapse of a massless scalar held. Perturbing around a 
continuously self-similar critical solution (the Roberts solution), we obtained an eigenvalue 
problem for the spectrum of perturbations. The remarkable feature of this model of critical 
scalar held collapse is that it allows an exact analytical treatment of the perturbations as 
well as of the critical solution, due to the highly symmetric background. 

An exact analysis of the perturbation eigenvalue problem reveals that there are no grow¬ 
ing non-spherical perturbation modes. However, there are growing spherical perturbation 
modes. Their spectrum is continuous and occupies a big chunk of the complex plane p 


In view of these hndings, the following picture of dynamics of scalar held evolution near 
self-similarity emerges: As we evolve generic initial data which is sufficiently close to the 
critical Roberts solution, non-spherical modes decay and the solution approaches the spher¬ 
ically symmetric one. Asymmetry of the initial data does not play a role in the collapse. 
The growing spherical modes, on the other hand, drive the solution farther away from the 
continuously self-similar one. In this sense, the critical Roberts solution is an intermediate 
attractor for non-spherical initial data. 

An interesting question, which is not answered by perturbative calculations, is the further 
fate of the scalar held evolution as it gets away from the Roberts solution. In all likelihood, 
it evolves towards the discretely self-similar Choptuik solution, which is a local attractor of 
lower codimension (one), as the continuous self-similarity is broken by oscillatory growing 
modes. After staying near the Choptuik solution for a while, the scalar held will eventually 
disperse or settle into a black hole, with these hnal states being global attractors in the phase 
space of held conhgurations. This evolution from attractor to attractor in phase space is 
somewhat analogues to a ball rolling down the stairs, going from a step to a lower step, until 
it reaches the bottom. 

The results of this paper shed some light at the complicated problem of critical collapse 
of generic initial data from the analytical viewpoint, conhrming the hypothesis that critical 
phenomena are not restricted to spherical symmetry. Investigation of the fate of a scalar 
held solution as it breaks away from continuous self-similarity, as outlined above, will further 
our understanding of the dynamics of scalar held collapse, and presents an interesting (and 
challenging) analytical problem. Numerical simulations might also help to establish a clearer 
picture of near-critical scalar held evolution. 
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APPENDIX A: PERTURBATIONS IN LONGITUDINAL GAUGE 


In this appendix we collect expressions for components of the perturbed Einstein-scalar 
equations calculated in longitudinal gauge (Jia = G = 0). The perturbed metric in longitu¬ 
dinal gauge is 

ds^ = huuYdu^ — 2(1 — huvY)dudv -f h^Ydv^ + (1 + KY)Pdil‘^, (Al) 


and the perturbed scalar field is 



+ FY. 


ui 


(A2) 


The Einstein equations for scalar field are equivalent to the vanishing of the tensor = 
R^i, — 20^^0 Its non-trivial components, calculated to the first order in the perturbation 
amplitude using the above metric and scalar field, are 


Ruu 


1 r 

2 - 


2('U Y ^hiiu,u 

— AvFy^ + /(/ -|- l)h^i 




uv 


(A3a) 


El,,, 


1 r 
'2 - 


("U UV^{huu,vv ‘^huv,vu T 

^huu,v “ 1 “ ( 2 "?/ H“ ^,v'} 

+ 2vF^y — 2uF^u — l{l + l)h^ ^ 


— uv 


(A3b) 


Rvv 


1 r 

2 - 


— 2( m ^ — uv)K^^^ + 2uhuv,v 

uh^'u^u (2u v^hyy^y Y 2uK^^ 
T Y 

Y 4:UFy Y Y l)hyy — -, 

J — uv 


(A3c) 


Ru9 2 UV^(^huu,v huv,u T ^,u) 


Y (2m — v)huv + 2vF 


U, 


e 


J u^ — uv 


(A3d) 


12 





+ {2u v)huy + 2vF , 

^ — uv 

(A3e) 

1 r 

FyO ~ (n ^)( ^UV,V T ^VV,U T 


1 Yff 

h 2F 

Ibyy , 

^ u — V 

(A3f) 

1 r 

Fy^p ~ (n ^)( hyy^y T hyy ^y T ^ 


Y 

h 2F 

IlUV ’) 

^ U — V 

(A3g) 

E/00 ~ 2(^1/- 


+ {2u - v){hyy^u + 2K^y) - 2uK^y 
- 2huv + 2/i^^ -2K + /(/ + l)i^] Y 


Y hyyY^QQ^ 

(A3h) 

= sin^ 9 Eoe, 

(A3i) 

EQ(P 9yyy(YQ^ COt 0 Y. 

(A3j) 


The scalar field equation requires the vanishing of □(;/), which, calculated to first order using 
the above metric and scalar field, is 

UV^F^yu T '^hyy^y 

+ uK^u — vK^y + 2uF^u — 2{2u — v)F^y 

- 2l{l + 1)f1 ■ (A4) 

J — uv 
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FIGURES 




(null coordinates) (scaling coordinates) 

FIG. 1. Global structure of the Roberts solution: The scalar field influx is turned on at u = 0; 
spacetime is flat before that. Field evolution occurs in the shaded region of the diagram, and there 
is a null singularity in the center of the spacetime. 
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